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ABSTRACT 

We present the results of a systematic search for gravitationally-lensed con- 
tinuum Lyman break 'drop-outs' beyond a redshift 7 conducted via very deep 
imaging through six foreground clusters undertaken with the Hubble and Spitzer 
Space Telescopes. The survey has yielded 10 z-band and 2 J-band drop-out can- 
didates to photometric limits of Jno —26.2 AB (5a). Taking into account the 
magnifications afforded by our clusters (1-4 magnitudes), we probe the presence 
of z >7 sources to unlensed limits of Jn —30 AB, fainter than those charted in 
the Hubble Ultradeep Field. To verify the fidelity of our candidates we conduct a 
number of tests for instrumental effects which would lead to spurious detections, 
and carefully evaluate the likelihood of foreground contamination by consider- 
ing photometric uncertainties in the drop-out signature, the upper limits from 
stacked IRAC data and the statistics of multiply-imaged sources. Overall, we 
conclude that we can expect about half of our sample of z-band drop-outs to 
be at high redshift. An ambitious infrared spectroscopic campaign undertaken 
with the NIRSPEC spectrograph at the WM Keck Observatory for seven of the 
most promising candidates failed to detect any Lyman a emission highlighting 
the challenge of making further progress in this field. While the volume density 
of high redshift sources will likely remain uncertain until more powerful facilities 
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are available, our data provides the first potentially interesting constraints on 
the UV luminosity function at z ~7.5 at intrinsically faint limits. We discuss 
the implications of our results in the context of the hypothesis that the bulk of 
the reionizing photons in the era 7 < z < 12 arise in low luminosity galaxies 
undetected by conventional surveys. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
formation — galaxies: high-redshift — gravitational lensing 



1. Introduction 



Very little is currently known about the abundance and luminosity distribution of star- 
forming sources beyond z > 7. The two principal techniques used to loc ate distant star 



form ing sources at lower redshifts, the Ly man-break 'drop-out' technique (Bouwens et al. 



2006 ) and the location of Lyman a emitters ( Kashikawa et al. 2006 : Shimasaku et al. 20061 ) 



become challenged by the lower performance of infrared instruments. In addition, the likely 
sources are much fainter, particularly if an incre asing fraction are sub-luminous as might be 
expec ted given mass assembly is at an early stage ( Loeb &: BarkanalboOl ; Choudhury fc Ferrara 
20071 ). Despite these hurdles, it seems reasonable to expect that there is an abundance of 
star-forming galaxies at these epochs. The improved measurement of the optical depth 
to electron s cattering derived from temperature-polarization correlations in the microwave 
background (IKomatsu et al.l 120081 ) suggests reionization occurred around z reion = 10.8 ± 1.4 
assuming it happened ins tantaneously; more probably it proceeded over an extended red- 
shift window 7 < z < 12 (ISpergel et al.l 120071 ) . Moreover, the detection of galaxies at z ~ 6 



with significant stellar masses and mature stellar populations (IStark et al.ll2007al : lEyles et al. 



20071 ) and the ubiq uity of ionized carbon in the intergalact ic medium probed by the high- 
est redshift QSOs (jSongailal 120041 ; iRyan- Weber et al.ll2006l ). together demand a significant 
amount of star-formation at earlier times, possibly enough to cause reionization. Although 
uncertainties remain, these independent arguments strongly motivate the search for z >7 
star-forming sources. 

Most of the early progress in this quest has been made through the publicly-available 
deep Hubble Space Telescope ( HST) images. The Great Observator ies Origins Deep Sur 



vey (GOODS. Dickinson et al.l (12003 )). the Hubble Ultra Deep Field (IBeckwith et all 12006 



UDF) and its associated parallel fields h ave b een used to search for 'drop-outs' in the 



(Bo uwens et al. 



2006 



Bunker et al.ll200 4). Zg50 ( IBouwens &: Illingworthll2006l : iBouwens et al. 



20081 ) and Jim ( Bouwens et al. 2005 . 2008 ) bands, corresponding to effective source redshifts 



of z~6, 7.5 and 10. These studies found a highly uncertain number density of candidates, 



- 3- 



none of which has been confirmed spectroscopically at z > 7. However, taken at face value, 
the overall conclusion from these ultradeep images is that the declining abundance of lumi- 
nous star forming sources beyond z ~7 is insufficient to account for reionization. Although 
there is no guarantee that star-forming sources did reionize the Universe at z ~10, a possible 
solution is that the bulk of the ear ly star formation r esides in an undetected population of 
intrinsically sub-luminous sources (IStark et al.ll2007al ). 



Prior to the availability of the next generation of telescopes, gravitational lensing is an 
effective means to evaluate this hypothesis. Depending on the method, foreground massive 
clusters can provide a magnification boost of x5 — 30 in fl ux (fo r unresolved sources) or 
in size (for resolved sources). An analysis by IStark et al.l (j2007d ) concluded that lensing 
surveys should be able to find ample candidates at z > 7, permitting follow-up spectroscopy 
and de tailed studies at s e nsitiv ity limits that would be unachievable otherwise. As pointed 
out by iBroadhurst et al.l (119951 ). this gain is offset by a reduction in the sky area surveyed, 
producing an overall increase or decrease in the number of lensed sources, depending on the 
slope of the luminosity function. 

A first atte mpt at constra i ning t he abundance of lensed drop-out galaxies at 6 < z < 
10 was made by iRichard et al.l (120061 ) using deep ISAAC near-infrared images obtained at 
the ESO VLT. A number of faint (intrinsic H(AB) ~ 26) candidates were identified in 
two clusters, where the magnification factor fi ranged from 0.4 to 2.5 magnitudes. This 
an alysis suggested a star formation rate density at z ~7.5 ~ xlO higher than that derived 



by iBouwens et al.l (120041 ) . Deeper NICMOS images have failed to confirm some of these 
sources. The number densities of faint candidates are currently being revised using follow- 
up imaging and spectroscopy for a small fraction of the ISAAC field of view (Richard et al., 
in preparation). The deeper NICMOS imaging undertaken in the present study illustrates 
the difficulties in making progress beyond z ~7 using ground-based facilities. 



In a parallel effort, IStark et al.l (l2007bl ) concentrated on the much smaller regions of very 
high magnification (// > 3 mags) termed the 'critical lines' of the lensing clusters. Using 
NIRSPEC on the Keck telescope they undertook a 'blind' spectroscopic search for lensed 
Lyman-a emitters in the redshift range 8.5 < z < 10.4. Despite the very small volumes 
probed in this unique survey, 6 faint candidates emerged across 9 clusters. Exhaustive follow- 
up imaging and spectroscopy has, so far, been unable to provide unambi guous confirmation 
of the nature of these sources. Taking into account the uncertainties, IStark et al.l 12007b 
concluded that the abundance of low-luminosity emitters in this redshift window may exceed 
0.2 Mpc -3 , suggestive of a major contribution of low luminosity star-forming galaxies to 
cosmic reionization. 



The caveats concerning this conclusion were discussed in detail by Stark et al. Although 
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representing a unique search for early star-forming sources at limits well beyond those probed 
otherwise, the volumes addressed are modest and significantly affected by cosmic variance. 
The 6 candidates were found in only 3 of the 9 lensing clusters; 6 clusters had no convincing 
candidates. Although Stark et al. were unable to prove, unambiguously, that the detected 
emission is due to Lyman a, the null detection of associated lines was used as an argument 
for rejecting lower redshift emission for most of their candidates.. Their conclusion that 
the bulk of the reionizing photons arose from low luminosity (~ 0.1 M yr _1 ) star-forming 
sources can be verified by this independent search for lensed continuum drop-outs at z >7. 

In searching for faint len sed drop-outs, t he ad vantage of HST over a ground-based survey 
such as that undertaken by iRichard et al.l (120061 ) is considerable. The ACS and NICMOS 
cameras are much more sensitive, not only because of the reduced background level relative 
to that produced by the atmosphere and by ambie nt temperature optics, but also because 
the typical sources have angular sizes of 0.2" or less (jEllis et al.ll200ll ). With similar exposure 
times, NICMOS can readily attain a depth of 26.5 AB mag, ~ 1 magnitude deeper than 
the earlier VLT/ISAAC project. Viewed through a typical z ~ 0.2 cluster, the NICMOS 
field (0.8 x 0.8 arcmin) closely matches the area of moderate to high magnification factors 
(yU ~ 1—4 mag). The effective increase in sensitivity to faint sources provided by the lensing 
magnificatio n along the line-of-sig ht to each NICMOS field probes limits fainter than those 
in the UDF (IBouwens et al.ll2008l ). albeit over a considerably reduced area. 



Early s t udies of lensed drop-outs with HST have served to illustrate the potential. 
Kneib et al.l (120041 ) located a source at z ~6.8 behind the cluster Abell 2218. This source 
forms a triply-imaged system with two bright elongated arcs, easily recognized as morphologica lly- 



simil ar in ACS and NICMOS images. Follow-up observations with Spitzer (jEgami et al. 



20051 ) provided improved constraints on the photometric redshift, stellar mass (M ~ 10 9 
M Q ) and past star formation history Very recently, a similar z ~7.6 lensed source was 
found by iBradley et al.l (120081 ) but this was not multiply-imaged. A semi-analytic analy- 
sis, empirically calibr ated using the lum inosity functions of Lyman-a emitters and drop-out 
galaxies at z ~ 5 — 6 (jStark et al.ll2007d ) predicts that such a NICMOS/ACS lensed imaging 
survey should typically detect 0.5-1 sources per cluster in the redshift range 7.0< z <8.5. 

The present program represents the logical next step: a concerted effort to verify the 
hypothesis advocated by Stark et al. as well as its associated predictions via deep imaging 
of 6 lensing clusters with HST (GO 10504: PI: Ellis) and Spitzer (GO-20439, PI: E.Egami). 
The primary goal is to determine the abundance of intrinsically sub-luminous z and Jno 
drop-outs, and to derive constraints on the possible contribution of low luminosity sources 
to cosmic reionization, independently of Lyman-a searches in blind or narrow-band surveys. 



The paper is organized as follows. In §2 we present the new HST/Spitzer and associated 
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ground-based imaging observations and their data reduction. A discussion of the means of 
selecting the drop-out candidates is presented alongside a catalog in §3. Issues of complete- 
ness and contamination by spurious and foreground sources are discussed. We conclude that 
a significant fraction of our candidates are possibly at high redshift. We then describe a 
Keck spectroscopic follow-up campaign in §4 which attempts to detect Lyman a emission 
from some of the most promising candidates. In §5 we discuss the UV luminosity function 
at z ~7.5 and review the implications in the context of the possible contribution of low 
luminosity to cosmic reionization. Our conclusions are summarized in §6. 

All m agnitudes given in this paper are standardized within the AB magnitude system 



(jOkd Il974h . We assume a flat universe and (O m , f^A, h)—(0.3, 0.7, 0.7) whenever necessary. 



2. Survey Strategy, Observations and Data Reduction 

We begin by discussing the selection of lensing clusters, the imaging datasets we have 
secured to select the various drop-out candidates and validate their high redshift nature, 
and the image processing steps taken to produce photometric catalogs. In addition to the 
HST and Spitzer datasets which form the fulcrum of this study, we have added ground-based 
imaging in both the optical and K band. In general terms, the HST data acts as the primary 
source of drop-out candidates and the associated data is used to constrain the likely redshift 
distribution. 



2.1. Lensing Cluster Sample 

Our primary criterion in selecting foreground clusters for this lensing survey is the value 
of the magnification factor expected for sources beyond z ~7 and the uncertainty implied 
in this estimate based on an understanding of the mass model. We considered a number of 
massive clusters at z ~ 0.1 — 0.5 with well-understood mass models capable of producing 
regions of strong magnification which match the imaging area of the HST cameras. Even 
though the area enclosed by the critical line for a z ~7 source is sufficient for the design of 
the program, a precise mass model is a clear advantage in determining the magnification of 
a particular source, as well as for predicting accurately the location of any multiple images. 

Our strategy parallels that discussed in some detail by Stark et al. Indeed that Lyman 
a critical line survey placed greater demands on the reliability of the cluster mass models 
as the magnifications are more extreme and positional uncertainties in the image plane are 
critical. Notwithstanding this challenge, Stark et al. found that the typical magnification 
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uncertainties for their candidates were only ~30% and that errors in the critical line location 
were usually only ±1 arcsec or so. 

In the present survey, six clusters was considered th e minimum number necessary to 
mitigate the effects of cosmic variance (IStark et al.ll2007d ) while being consistent with the 
observing time available. In considering the final tally of clusters, we included clus ters which 
we have modeled with the latest version of the Lenstoofj] software f Julio et al. 2007 ). This im- 
proved code provides a new Bayesian optimization method to der i ve err or es timates for each 
optim ized parameter. Following the discussion in iRichard et al.l (120071 ) and iLimousin et al. 
( 120071 ). this optimization can be used to compute error estimates on the individual magnifi- 
cation factors for each of our sources. 

Wherever possible, we included clusters with usefully deep archival HST and Spitzer/IRAC 
data. Deep optical (AB>27.0) ACS and/or WFPC2 images, previously used for the iden- 
tification of multiple images during the development of the mass models, allow low redshift 
contaminants to be identified reliably. IRAC images at 3.6 and 4.5 fim, to the same 1 /xJy 
sensitivity as our previous work in Abell 2218 lEgami et al.l (120051 ). are available as part of a 
Spitzer Lensing Cluster Survey (GTO 83, PI: G. Rieke) or publicly from the archive (GTO 
64, PI: G.Fazio). Finally, we required that each of the selected clusters be visible from 
the northern hemisphere, in order to facilitate spectroscopic follow-up and further i^-band 
imaging with the Keck and Subaru telescopes. 

The six clusters satisfying the above criteria are presented in Table [TJ Although the 
references cite the most recently-published mass models, as described above, in each case 
we have utilized the available multiple images and their redshifts in improving these mass 
models using Lenstool. Four out of the six clusters are in common with the sample adopted 
by Stark et al. 



2.2. Hubble Space Telescope Data and Reduction 

Our large program with HST (GO 10504, PI: Ellis) comprised deep observations in the 
-2850, ^ioo and Hieo bands, using the Wide Field Camera of ACS and the NIC3 configuration 
of NICMOS. The region enclosed by the critical line for the putative z >7 sources was 
typically covered by two NICMOS pointings per cluster, usually adjacent. This ensured a 
magnification gain of \i ~ 1 to 4 magnitudes, with a typical value \i ~ 2 mags, throughout 
the NICMOS imaged field (Fig. [[]). The total sky area covered by the NICMOS observations 



1 Publically available, see 



http : //www . oamp . f r/ cosmology/lenstool 



to download the latest version 



- 7- 



is 8.9 arcmin 2 for the six clusters. 



The ACS images were reduced with the multidrizzle software (jKoekemoer et al. 



20021 ). This removes cosmic rays and bad pixels and combines the dithered frames to correct 
for camera distortion. The output pixel scale was fixed at 0.04"and we used a pixf rac value 
of 0.8 for reducing the area of the input pixels. We made small corrections to the absolute 
astrometry to allow for ACS frames taken at different epochs (e.g. Abell 2218 in the F850LP 
band, see Table [2]). These corrections were computed by correlating the catalog positions of 
bright objects detected in the overlapping regions. 

Each single set of NICMOS observations consists of 8 (in F110W) and 10 (F160W) 
frames of ~ 1000 sees, taken with the NIC3 camera using the SPARS64 or SPARS 128 
sampling sequences. A basic reduction was performed by adopting the procedures given in 
the NICMOS data reduction handbook^. Starting with the post-calibrated frames, bad pixels 
are flagged and rejected based on individual histograms, cosmic ray are rejected using the 



LAC0SMIC (Ivan Dokkuml 120011 ) IRAF procedures, frame-to-frame shifts are measured using 
a cross correlation technique, and all frames are drizzled onto a NICMOS pixel scale (0.2") 
to produce an initial reduced image. This then serves as a comparison for each individual 
frame so that deviant pixels can be flagged to improve image quality in a second drizzling 
operation, this time undertaken with an image pixel of 0.1", to obtain a better sampling. 

These initial images reveal a number of cosmetic effects (bias and flat residuals, bad 
columns and bad pixels, quadrant-to-quadrant variations, background variations) that led 
to a second stage of image reduction. We used an improved pixel mask to flag several bad 
columns and bad pixels close to the frame edges, and we examined each frame individually 
in order to remove bias and flat residuals and quadrant-to-quadrant variations. Finally, we 
subtracted a smoothed background obtained by averaging the frames of all the observed 
clusters, masking every pixel lying on a physical object. In order to combine all NICMOS 
observations of a given cluster, usually taken at different epochs and with slight variations 
in the sky orientation, the individually-reduced images were registered onto the wider field 
ACS data prior to the final drizzling procedure and combination into a NICMOS mosaiced 
image. 

The NICMOS data acts as the primary basis for selecting our drop-out candidates in 
association with non-detections in very deep optical data. In addition to our own ACS data, 
we reduced deep ACS and WFPC2 images from the archive for each cluster in our sample 
(Table [2]). As with the ACS F850LP data undertaken in our own program, we reduced these 
data using the IRAF procedures multidrizzle and drizzle as discussed above. 



^http : //www. stsci . edu/hst/nicmos/documents/haiidbooks/haiidbooks/DataHandbookv7/ 
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2.3. Ground-based Data and Reduction 

As mentioned, optical non-detection to deep limits is a key necessity in considering the 
validity of our drop-out candidates. Most clusters have complete optical coverage with HST 
(from ~ 400 to ~ 850nm ) useful for this purpose. However, Abell 2219 has F702W only, and 
Abell 2390 has F555W and F814W only. For these two clusters, additional deep (R ~ 26.7 
and / ~ 25.6 at 5a in 1.2"diameter aper tures) ground-based images taken with the CFH12k 



camera on CFHT (IBardeau et al.l 120071 ) were also examined to check for non-detections at 



the locations of the drop-outs. 

Ground-based photometry in the K band can likewise provide additional information 
for z drop-outs, improving the photometric redshift estimates, and for the reliability of J 
drop-outs where otherwise only a single color would be available (see Section I37L]) . Although 
a challenging undertaking given the depth of the ACS and NICMOS data, we conducted 
various f^-band imaging observations of clusters using the Keck and Subaru telescopes. 

The Near Infrared Camera (NIRC) on Keck-I was used on July 2006 to observe the 
central regions of the clusters Abell 2218 and Abell 2219. The square field of NIRC (38"a 
side) is slightly smaller than that of NICMOS therefore we concentrated these observations 
on our best z and J-band drop-outs in each cluster (see section l3TTj) . The seeing was stable 
and in the range 0.5-0.6". We used dithered exposures of 6 coadds x 10 sees exposure time 
each. The NIRC pointing in Abell 2218 was partially covered, in the region of one of the 
drop-outs, by previous NIRC observations obtained in 22/23 July 2002 (A. Blain & N. Reddy 
2003, private communication). 



The Multi-Object InfraRed Camera and Spectrograph (MOIRCS) (jlchikawa et al.ll2006l ) 
at Subaru was also used during two observing campaigns, in August 2006 and May 2007. The 
larger field of view (7' x 4') of this camera ensures a complete coverage of the NICMOS, ACS 
and IRAC data in each cluster. The four clusters in common with the Stark et al. sample 
(Abell 2218, 2219, 2390 and CL0024) were imaged for ~ 5 hours each under very good seeing 
conditions (0.3-0.4") using dithered exposures of 50 sees duration. We used the MCSRED 
software package^ to perform the flat-fielding, sky-subtraction, distortion correction and 
mosaicing of individual images. The MCSRED package also includes MOIRCS-specific tasks 
which correct for quadrant shifts and to fit sky residuals. The depth reached by the MOIRCS 
imaging data [Kab ~ 26.1) in Abell 2218 and ABell 2219 is similar or deeper than the NIRC 
observations of the same fields. Therefore, we used the NIRC images as independent check 



3 available from |http://www.naoj.org/staff/ichi/MCSRED/mcsred. html 
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for the MOIRCS photometric measurements, performed under better seeing conditions. 



2.4. Spitzer Data and Reduction 

Each cluster has been observe d at 3.6, 4.5, 5.8, and 8.0 fim using the I nfrared Array 



Cam era (IRAC; iFazio et al.l (120041 )) on board the Spitzer Space Telescope ( jWerner et al. 
20041 ) . In this paper, we only discuss the first two channels (3.6 and 4.5 fim) where the 
depth achieved is potentially useful in comparison with our ACS and NICMOS data. Each 
of the IRAC channels uses a separate detector array, and the 3.6 /zm channel (A c = 3.56/xm; 
AA = 0.75/im) and 4.5 /im channel (A c = 4.52/im; AA = 1.01/zm) use 256 x 256 InSb arrays 
with a pixel scale of 1'.'2 pixel -1 , producing a field of view of 5'.2x5'.2. A frame time of 
200 seconds was used with the small-step cycling dither pattern initially, but this was later 
changed to 100 seconds with the medium-step cycling pattern for a better removal of cosmic 
rays and other artifacts. 

Most of the IRAC data come from the Spitzer GTO program PID:83 (PI: G. Rieke). 
The total integration time is 2400 s per channel, usually obtained at two different epochs. 
CL0024 was observed as part of another GTO program PID:64 (PI: G. Fazio) with the 
integration times of 2400 s at 3.6 /j,m and 3600 s at 4.5 /im. Abell 2218 was the subject of a 
deeper GO campaign (PID:20439, PI: Egami), and was observed for ~10 hours per channel, 
split into 6 separate observations (AORs). When these GO data are combined with those 
of the GTO program (PID:83), the resultant total integration time is 37700 seconds per 
channel. The corresponding depth (25.5 AB) in this particular cluster is much closer to the 
magnitude of our high redshift candidates, (see Table [2] for a summary of the data). 

In producing the IRAC images, we started with the basic calibrated data (BCD) of 
each individual frame produced by t he SSC pipeline, an d combined them using a custom 



IDL mosaicing routine as presented in lEgami et al.l (120051 ) . The final pixel size of the IRAC 



mosaics is 0.6 "/pixel, half that of the original data. A conservative estimate for the absolute 
calibration uncertainty is 10%. 



2.5. Foreground Subtraction 

Although the magnification afforded by lensing clusters offers a unique gain in prob- 
ing the distant universe, the central regions are dominated by bright, extended spheroidal 
galaxies which obstruct and whose light increases the background level. More importantly, 
it also affects the photometry and color measurements of any underlying fainter source. For- 
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tunately, the majority of these galaxies are also good light deflectors, their morphology is 
usually regular and their light distribution can be accurately modeled by a sequence of ellip- 
tical isophotes. Since the major contributor to this deficiency is the brightest cluster galaxy 
(BCG), we have modeled and subtracted its light for each cluster in the ACS (J 8 5o band) 
and NICMOS images using the IRAF task ellipse. Apart from a small region (~ 1.0" 
radius) around the core of the removed BCG, this frees the image from the majority of the 
contamination in highly magnified regions close to the critical line, and aids in the detection 
of underlying sources (see Figure [2]). 

Foreground subtraction is more challenging for the IRAC images because of the much 
coarser PSF. The smaller number of available pixels per galaxy prevents the use of the 
ellipse modeling technique discussed above. Instead we adopted the ellipse model of 
the BCG derived from the NICMOS data and, on the assumption that the spectral energy 
distribution is identical across the galaxy at these long wavelengths, convolved this model 
with the IRAC PSF derived by stacking ~ 50 bright isolated point sources. This model 
of the BCG was subtracted after an appropriate scaling factor. We found this procedure 
to be very effective, with residuals from the subtraction confined within a 2.0"radius region 
(Figure [2]). 



2.6. Final Photometry 

The primary filters we will use for identifying the near-infrared drop-outs are the ACS 
z 850 and NICMOS Jn , Hi 60 bands available for each cluster. The bulk of the rest-frame UV 
flux is contained in the two reddest filters at high redshift, so we require all sources to be 
detected in both. The signal to noise ratio is improved by combining the J and H images, 
once normalized to a similar noise level, to form a single < J + H > detection image. This 
provides a more accurate measurement of the centroid and geometrical parameters of these 
objects, 



We use the "double-image" mode of the SExtractor package iBertin fc Arnoutd (119961 ) 
to detect objects and compute magnitudes within a 0.6"(NICMOS images) or 0.3"(ACS 
images) diameter aperture. Corrections to total magnitudes, assuming a point source, were 
estimated using bright isolated unsaturated stars. The corresponding values are 0.3 and 0.6 
magnitudes for ACS and NICMOS, respectively. 



The drizzling procedure used in the HST reductions, while conserving flux, does intro- 
duce correlations between neighboring pixels and hence unreliable error estimates (jCasertano 
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2000l ). We applied their equation [A20] to the SExtractor photometric errors. Dithered expo- 
sures also introduce a varying effective exposure time across the mosaic, this effect is clearly 
apparent close to the boundaries of the field. We constructed a corresponding weight map, 
proportional to the effective integration time at each pixel position in the detection image, 
and used it as an input parameter of SExtractor. This ensures source detection at a uniform 
noise level across each image. 

The background noise level a was measured for the ACS and NICMOS images in order 
to estimate the achievable photometric limits in each band. The 5a values are reported in 
Table [2] as total magnitudes. We computed the completeness limits in each NICMOS filter 
used for source detection by adding artificial unresolved sources in the magnitude range 23-28 
AB. Such sources were added 1000 times at 30 different random locations on the image and 
then extracted using the same photometric techniques as applied bf to the science image. In 
this manner we compute the point source completeness as a function of source magnitude. 
Only blank regions were chosen for this exercise, defined after application of a 5cr threshold 
to mask pixels associated with bright objects. The 50% completeness values derived in this 
way are listed for each cluster and band in Table El Typically, our NICMOS data is 50% 
complete to Jn ~25.9 AB and to H 160 ~26.05 AB. For each eventual drop-out candidate 
(§3.3), we will assign a completeness factor, S comp , based on its magnitude, that will be taken 
into account in estimating the high redshift luminosity function. 

Our NICMOS survey reaches limits of typically 26.2 and 26.5 AB in J and H, re- 
spectively, over 8.9 a rcmin 2 . For comparison, the UDF limits are 27.8 AB over 7 arcmin 2 



(IBouwens et al.ll2004l ). However the magnification provided by the foreground clusters (Fig. 
[I]) enables us to reach unlensed source magnitudes (assuming z ~ 7) of 28-30 AB over a total 
area of 0.1-1 arcmin 2 . 



3. High Redshift Candidates 

We now turn to the selection of our high redshift candidates. The primary selection will 
be based upon the zsso — Jno color for the z drop-outs at redshifts z ~7-8, and the Jho — #160 
color for the J-band drop-outs at redshifts z ~8-10. This section has two components. First 
we discuss the optimum color criteria for drop-out selection, the degree of completeness and 
issues of possible contamination from spurious artefacts at these faint limits. In this way 
we establish a robust set of candidates whose likelihood of being at high redshift we then 
assess in the second part of this section using additional criteria including their photometry 
at other wavelengths and location with respect to the cluster mass model. The catalog of 
candidates is summarized in 53.3. 
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3.1. Photometric Selection and Completeness 

The primary concern in selecting high redshift drop-outs from photometric data alone 
is the issue of contamination from lower redshift objects, including z ~2 early-type galaxies, 
dust-reddened objects over a wide redshift range or low mass Galactic stars with deep molec- 
ular absorption bands. Figure [3] illustrates the problem. A single color-cut fails to isolate 
z-band drop-outs from a variety of z ~2-4 sources and the confusion is worse for J-band 
drop-outs (Fig. @J. 

This problem has formed the basis of much discussion in the literatu re. For the z- 



band drop-outs, contamination can be reduced by considering a secon d color flStanway et al. 



20041 ). of which the most useful with NICMOS data is J no - H im flBouwens et all 120041 ) . 
Star- forming galaxies at high z should display a prominent discontinuity in z — J while 
remaining blue in J — H, as illustrated in Fig. [3j The first issue we address is the optimum 
cut in both colors, on which depend both the redshift range explored and the amount of 
contamination by lower redshift objects. 

Figure [3] shows color-redshift tracks for various galaxy classes (and also includes, for 
convenience, the colors of our eventual candidates discussed in §3.3). The location of these 
tracks suggests the following prescription for selecting sources with redshifts 6.8 < z < 8.0: 



(^850 - Jno) > 0.8; (z 850 - J 1W ) > 0.66 (J no - H lm ) + 0.8; (J no - H im 



< 1.2 



The baseline (z — J) > 0.8 color selection we adopt above is similar to that used by 
Bouwens et al.l (120041 ) and lBouwens et al.l (120081 ) . However, as the photometric errors for our 
candidates are typically 0.2-0.3 magnitudes, the probability of low-redshift contamination 
remains significant in the range 0.8 < (z — J) < 1.1. We also ex plore the use of a more 
restric tive color cut (z — J) > 1.25, close to the criterion used by iBouwens fc Illingworth 
(120061 ) ((z — J) > 1.3), also shown on Figure [3j 



Application of this color selection reveals 10 possible candidates of which only two 
satisfy the more rigorous (z — J)ab >l-25 color cut. We discuss the merits of each of these 
candidates in more detail in the following sections. A striking feature is that all have J-H 
colors bluer than predicted by the redshifting of local spectral energy distributions. Similar 
claims for blue restframe colors have been made for z-band drop-outs (Stanway et al 2004). 

Concerning the J drop-outs, corresponding to the redshift rang e 8.0 < g 5, 10- 5 , cuts 
of ( J — H) > 1.8 and (J — H) > 1.3 were adopted respectively by IBouwens et al.l (120051 ) 
and lBouwens et al.l (120081 ). In their shallower, larger area survey aimed at locating luminous 
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J drop-out candidates, iHenry et al.l (120071 ) adopted a more restrictive (J — H) > 2.5 cut. 
The large variation in these color cuts reflects the differing depths of optical exclusion in 
the various samples. The sample selected by Henry et al. is limited to NICMOS-detected 
sources in 2 bands only, with no deep supporting optical observations. Therefore, despite the 
apparent stringent color cut in (J — H), it is still more likely to suffer from contamination 
by lower redshift sources. Noting our deep optical data, we adopted a (J — H) > 1.3 cut, 
revealing two candidates (Fig. HJ). Neither would satisfy the (J — H) >1.8 cut. 

We now turn to the important question of how complete our likely drop-out sample will 
be given our adopted magnitude limits with ACS and NICMOS. We can easily imagine that 
genuine drop-outs will be missed because photometric errors will scatter points outside our 
selection region; likewise, lower redshift sources will be scattered into our color selection box. 

Using the procedures adopted to determine the 50% magnitude completeness limits in 
Section 12.61 we can estimate the color selection completeness and color selection contami- 
nation by introducing artificial objects with a range of J magnitudes (AB=24 through 27) 
whose z — J breaks have a range (0.6 through 1.2) and determining what fraction of objects 
which would lie outside of the box with perfect photometry but are scattered in, and what 
fraction of objects lying truly inside are scattered out. We assume a flat /„ spectrum for the 
UV continuum between the J and H bands so that, with reference to Figure 3, the problem 
becomes effectively one dimensional. 

The two panels in Fig. [5] show the results of this test. The selection completeness 
fraction f comp (left) represents the fraction of objects of a given magnitude and break color 
that we are able to recover in our selection. The contamination fraction f cont (right) accounts 
for the fraction of objects with a lower break (z — J< 0.8) that are photometrically scattered 
in the color-color diagram so that their observed magnitude and color would allow them to 
enter our selection window. As with our magnitude completeness function, S comp (§2.6), 
both of these correction factors will be used, for each candidate, to correctly estimate the 
true number density of objects having colors corresponding to z > 6.8 galaxies (i.e. a break 
z — J above 0.8). Individual correction factors are reported in Table [31 

Examining Fig. [5J as expected we find negligible difficulties for the brightest sources, 
but contamination and incompleteness become more troublesome at fainter magnitudes, 
depending on the (z — J) color. We find that the selection completeness ranges from 50 to 
95% and the contamination fraction is typically 15%. 
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3.2. Verifying the Candidates 

The selection techniques discussed above yields a list of candidates for more careful 
examination. Here we discuss further tests to determine the possibility that some might be 
spurious prior to establishing a catalog of genuine sources whose redshift distribution we will 
explore using our additional photometric data. 

3.2.1. Spurious Detections 

As the signature of both our z and J-band drop-outs consists of a non-detection in 
the optical band, we must seriously consider the possibility of spurious detections in the 
NICMOS data. This is particularly the case for the J drop-outs where only a single band is 
involved. An optical non-detection is defined as a implied flux lower than the photometric 
limit (5 a in a 0.3"-diameter aperture, for ACS and WFPC2 data, see Table [2]). These 
measurements were made using SExtractor in its "double-image" mode after the data was 
resampled and aligned to the NICMOS images. Further measures were made using the 
original multidrizzled images (see Sect. 13.51) . 

A visual inspection was performed for each candidate in order to reject obviously spu- 
rious detections in the NICMOS images or false non-detections in the optical bands. The 
astrometric position of the candidates was used to perform this examination on the original 
images, to prevent biases arising from resampling and geometrical transformations. During 
this process, we rejected a number of candidates due to their proximity to the center of the 
removed BCG, or due to obvious contamination from very bright galaxy haloes, both leading 
to noisier or biased photometric measurements. 

Our photometric detection is based on the combination of 10 individual NICMOS im- 
ages per pointing (4 and 6 exposures in the Jn and H 160 band, respectively). Because of a 
significant number of remaining hot and dead pixels in these individual frames, we investi- 
gated the fraction of spurious sources that would contaminate our photometric catalogs. To 
quantify this problem, we constructed a noise image for each cluster and near-infrared band 
whose purpose is to remove signals from all genuine sources while maintaining the same noise 
properties as the real data. This was done by subtracting in pairs the individual frames prior 
to eventual coaddition. We then applied our usual photometric detection software, using the 
same parameters as in the original images. 

This noise image reveals residuals near the frame edges (due to the dithering process 
undertaken during the observations) and in the centers of the brightest objects, which were 
masked out in the detection process (Figure [6]). The affected area accounts, in total, for 
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~ 20% of the NICMOS field. By comparing the number of spurious sources detected in 
the noise image with the number of objects present in the original catalogs, we estimate 
the spurious fraction in the magnitude range of our candidates (24.7 < J < 26.3 and 
25.3 < H < 26.7) to vary from 4% to 18% from cluster to cluster, with an average value 
of 10%. Typically, therefore, we can expect around 90% of our candidates to be robust 
astronomical detections. 



3.2.2. Detector Remnance 

One specific worry, not addressed in the tests mentioned above, concerns the possibility 
of image persistence or ghosts appearin g in the NICMOS fram e during, or shortly after, 



an overexposure by a very bright source (IDaou fc Skinner! 1 199 71 ). The source producing the 



largest number of counts (~ 25000 ADU) in our dataset is a Jab = 17.9 star in Abell 2390. 
We do not detect any measurable electronic ghost for this particular source, but persistence 
is seen at the level of a Jab — 24.9 spurious source in an exposure taken 20 minutes later. 
Once averaged over the entire sequence (6 exposures), this persistence corresponds to a 
Jab = 26.8 source, i.e. fainter than our detection limit. We also verified, for each pair 
of successive exposures, that none of our candidates is coincidently located at the relevant 
position of similarly bright (J < 19) stars in all clusters. 

Although NICMOS exposures from independent orbits are usually separated by a > 30 
minute delay, persistence might arise as a result of much brighter sources observed immedi- 
ately prior to execution of our program. To eliminate this possibility, we located all preceding 
exposures and found no measurable persistence in our data, even in the case of a calibration 
program (GO 10726, PI: de Jong) aiming at measuring non-linearity effects by repeatedly 
saturating the NIC3 detector. 



3.3. Catalog of Drop-out Candidates 

As a result of our visual inspection of the initial candidate list selected according to the 
precepts of Section 3.1, we emerge with 10 candidate z drop-outs and 2 J drop-outs. As seen 
in Fig. [31 2 of the 10 z-dropouts satisfy the most restrictive color-selection (z— J) > 1.25. The 
photometric measurements for these are summarized in Table [3] and the relevant detection 
images are presented in Fig. [71 Following the tests described above we can expect over 90% 
of these to be genuine astronomical sources. 
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3.4. Redshift Estimation 

We now turn to the important problem of contamination by lower redshift sources. We 
first use the spectral energy distribution (SED) of each candidate to estimate the individual 
photometric redshift. We then consider statistical arguments that can be applied to our 
entire candidate population. 



3.4-1- Photometric Redshifts 



As we have seen, the two color selection presented in Sect. 13.11 enables us to select 
high redshift galaxies with some confidence. However, we can use the multiwavelength data 
available for each source, including all upper limits arising from non-detections, in order to 
derive a photometric redshift probability distribution. To a ccomplish this, we use d an up- 
dated version of the photometric redshift software HyperZ (IBolzonella et al.ll2000l ). Best-fit 
redshift distributions were comp uted using a standard SED procedure wit h a variety of tem- 
plates, including e mpirical data (|Kinnev et al.lll996l ; IColeman et al.lll980l ) and evolutionary 
synthesis models ( jBruzual fc Charlotl 120031 ). We searched the redshift range < z < 10, 
whi le reddening w a s kept as a free parameter ranging between Ay = and 2 mag, assuming 
the ICalzetti et al.1 (120001) law. Th e effect of Lyman forest blanketting is included following 
the prescriptions of iMadaul (119951 ). 



This approach is only practical for the z drop-out candidates and the main results are 
presented in Fig. [7] (rightmost panels) where we overplot the best fit templates on each SED 
and present the redshift probability distribution P(z), marginalized over the e ntire parameter 
space (templates and red dening). As has been found by many workers (ILanzetta et al. 



19961 : iRichard et al.ll2006l ). the likelihood function reveals two peaks with different relative 
intensities, the relative height of the lower redshift peak acting as a valuable measure of 
foreground contamination, as it is linked to the color degeneracy between high-redshift star- 
forming galaxies and lower redshift early- type or reddened objects. 

Encouragingly, each candidate is more likely to be at high redshift and the probability 
of a foreground source is negligible for 4 out of our 10 sources. Integrated over all our 
candidates, we use P(z) to compute the probability that each object lies beyond a redshift 
of 6, «6 = P{z > 6.0), after normalizing P(z) to unity over the redshift range < z < 10. 
We find «6 values ranging from 0.46 to 0.91, with an average value of 0.65. 
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3.4-2. Stacked Imaging 

As is evident from Table [31 the individual candidate SEDs are mostly restricted to 
detections in only one or two bands, with upper limits at other wavelengths. Although this 
precludes precise photometric redshift measures for each candidate, we can make further 
progress by combining data over several adjacent wavebands, and also by examining the 
average SED of the population to see if it is statistically consistent with that expected for a 
high redshift source. 

Several of our clusters have ACS data in multiple bands (Abell 2218, CL1358, C10024) 
which we combined after aligning the images with integer ACS pixel shifts and normalization 
each to a constant signal to noise ratio. The depth of this combined image, covering the 
wavelength range 4500 - 8000 A , is typically 0.6 to 1.0 magnitudes deeper than the individual 
bands. Yet in each of the 5 z drop-outs, no optical detection is seen, increasing our confidence 
level in the corresponding drop-outs. 

We likewise generated stacked z, J and H images for all 10 z drop-outs, as well as a 
stacked K image for the 4 sources observed with MOIRCS. In this case, We selected a 10" x 
10" region around each candidate, and averaged the data rejecting 20% of the outlier pixels. 
For 6 objects free from contamination from nearby objects, the same stacking procedure was 
performed in the 3.6 and 4.5 /zm IRAC bands. The stacked images are presented in Fig. [SJ 
No significant flux was detected in either the stacked MOIRCS or IRAC images. 

Using the SExtractor parameters adopted for processing the individual images, an av- 
erage SED was constructed for 3 populations depending on the availability of multi-band 
ACS, K band and IRAC imaging. The photometric properties of each are listed in Table H] 
and the SEDs are shown in Fig. [9j These various stacked datasets offers a new opportunity 
to address the question of foreground contamination. 

Using HyperZ, we fit these three SEDs exactly as described in Sect. 13.4.11 The optimal 
photometric redshift is consistent in each case (Fig. fTDl) . with z ph = 7.35 ± 0.07, and a 
reduced x 2 lower than 1. We also fitted the same SEDs, but now restricting the photometric 
redshift to the range < z p h < 3. The best fit at low redshift is found at z p h = 1.75 ± 0.03, 
but with a much higher x 2 , between 5 and 10. To the extent that the low redshift solutions 
are credible, they imply galaxies with typically ages of 500-700 Myrs and extinction values 
of A y ~ 1.0-1.4. 

One question that arises is whether the near-infrared J — H color for our z drop-outs 
is consis tent with expectation s, and also with that of the other limited detections at high 



redshift (IBouwens et al.ll2008l ). Using the stacked images, we find this color is typically very 



blue: Jiiq — Him ~ 0.0. We define the rest-frame ultraviolet slope (3 as f\ oc A 13 between the 
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restframe wavelengths 1500 and 2000 A ( Calzetti et al.lll994l ). and estimate the uncertainty 
in (3 for our sources either from the dispersion within the range of plausible redshifts (at ±1 
a), or by using the adopted photometric error bars in J and H. In both cases, we find a 
consistent value: (3 ~ 2.8^Q2 5 . The mean slope is somewhat blue r than the average val ue of 
2.0±0.5 for a sample of i band drop- outs measured at z ~ 6 by iBouwens et al.l (120061 ) but 
within the range of (3 ~ 3.0 found by lBouwens et al.l (12004 ) in their sample of z drop-outs, 
and m arginally consistent with (3 ^ 2.5 given for several objects by lBouwens &: Illingworth 
(120061 ). In order to produce such a blue ultraviolet slope, the SED can only be fit by very 
young models (typically < 100 Myrs) with little or no reddening (Ay < 0.1). 



3.4-3. Galactic Stars 



Our next test for for eground contamination concerns the questi on of low mass Galactic 



stars. Various authors ( Stanway et al.l 12004 ; IBouwens et al. 



20041 ) have pointed out the 



difficulty of using optical-near infrared colors to d istinguish betwee n cool stars and breaks 
arising from Lyman absorption at high redshift. L (ICruz et al.ll20071 ) and T (lBurgasserll2004l ) 
dwarfs exhibit metal and H2O absorption features which produce features simila r to the drop- 



out s ignature. Indeed, L and T dwarfs from the SDSS and 2MASS surveys (IKnapp et al. 



2004 ) lie well within our color-color selection region (Fig. [3]). Likewise, we find \ 2 values 
similar to the best template spectrum from HyperZ when fitting the SEDs of the candidates 
with a library of L0 and T dwarf templates spectra. 

In such circumstances, HST data has occasionally been used to evalua te whether the 



sourc es have half-light radii Rh consistent with being extragalactic objects (jStanway et al. 



20041 ). However, this is not a definitive criterion as we already kn ow that many spe ctroscopically- 
confirmed low luminosity high redshift sources are unresolved ( lEUis et al.l 1200 ll ). We mea- 
sured Rh, defined as the radius enclosing half of the flux in the detection (NICMOS) images, 
and compared the values with that derived for bright not-saturated stars (0.2"). Figure [TT1 
represents the location of our candidates in a Rh vs Jab diagram, together with all other 
objects in our photometric catalog. Unfortunately, surface brightness dimming generates 
a cut-off at large values of Rh, which does not allow to distinguish resolved and resolved 
sources at the faintest magnitudes. This is the case for 4 of the z-dropouts. At most we can 
say that 2 dropouts are unresolved, whereas 4 are resolved. 

Noting the difficulty of separating stars from galaxies at these faint limits, a more 
pr actical app r oach is to examine the likely contamination statistically. Using simulations 
by iBurgasserl (120041 ). we computed the number of expected low-mass stars in the total area 
surveyed with NICMOS (7.7 arcmin 2 ), for all spectral types between L0 and T8, up to the 
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extent of the thick disk (1 k pc). Assuming a slop e a = 0.0 for the mass function, consistent 
with recent observations by iMetchev et all (120081 ). we predict only 1.1 dwarf in our survey. 
A more pessimistic a = 0.5 slope yields 1.5 stars. Thus, while there is undoubtedly some 
uncertainty surrounding these predictions, it does seem unlikely that cool Galactic stars 
represent a significant contaminant at such faint limits. 



3.5. Magnification and Multiple Images 

Our final test concerning the high redshift nature of our candidate drop-outs concerns 
their location within the image plane of the lensing cluster. A key question is whether any 
mi ght be expected to b e multiply-imaged as was the case, for example, for the object studied 



by iKneib et all (12004 ) . Most of the area covered by our NICMOS observations lies within 



the region of high magnification (/i > 2.5 mags) where multiple images may occur. 

Using the mass models for each cluster (updated using Lenstool from those referenced in 
Tabled]), we estimate the magnification of each drop-out based on its photometric redshift 
and location and examine the likelihood of any counter images as well as their predicted 
location and relative fluxes. The Bayesian optimization method incorporated in Lenstool 
also provides the uncertainties in these quantities. 

In the majority of the cases (7 out of the 10 z drop-outs), the model predicts a pair 
of images with similar fluxes (within 0.2-0.4 magnitudes) straddling the critical line (Figure 
[T2l . Single images are expected in two other cases (CL1358z3 and A2667zl) and in the final 
case (CL1358zl), the objects sits on the predicted critical line, but is not expected to be 
viewed as a distinct pair at the finite angular resolution of NICMOS. These three cases are 
consistent with our observations. Considering the two J drop-outs, A2219jl is predicted to 
be another example of close merging unresolved by HST, and A2667jl is predicted to have 
much fainter counter-images (by 0.8 to 1.5 mags), below our detection limits. 

Our attention is thus focused on the 7 cases where second images are expected. The 
typical positional errors are around 1-3 arcsecs. Unfortunately, most are located outside 
the area surveyed by NICMOS or in regions close to the edge of the detector or under the 
central BCG where the noise level is high. Only in two cases, A2218zl and A2667z2, does the 
mass model predicts a detectable counter-image in a relatively clean region of the NICMOS 
detector. Unfortunately, no significant flux, within the range expected, was seen at either 
position in the J-band of H-band image. 



Unfortunately, this test, valuable in principle, provides an inconclusive outcome in our 
exhaustive quest to confirm the high redshift nature of our candidates. Although the high 
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redshift test is only applicable to two of our 10 z-band drop-outs, we are unable to see either 
of the counter images. However, an identical search for multiple images was performed 
under the assumption that each candidate is at a lower redshift, z ~2. We would then 
expect counter-images at different locations, closer to the cluster center. However, we did 
not detect any of these predicted counter images either, for 4 clear cases. Only in a more 
marginal J-drop candidate, which we dropped at an early stage, did this test succeed in 
demonstrating a z ~2 solution (Fig. 12). Given the strong likelihood that most of the 
candidates lie either at z ~2 or z ~7-8, the test has a confused outcome. On the one hand, 
for those high z candidates where we could expect to see a counter-image, none is seen. On 
the other hand for a larger sample of candidates, assumed to be at z ~2, none is seen either. 
Accordingly, we deduce the test is not effective as a redshift discriminant. 



3.6. Summary 

We now summarize the possible success rate of our survey in generating high redshift 
sources, concentrating on the z-band drop-outs. Out of our ten sources, we find that at 
most one source is spurious and ~one is a Galactic star. Thus we conclude 8 are likely to 
be extragalactic sources. Admittedly, about half of our candidates are unresolved, but we 
believe this is to be expected given the intrinsically faint limits we are probing with our 
lensing method. Discounting the inconclusive test based on counter-images (Section 3.5), 
and noting the 35% contamination from z ~2 sources, we conclude that, statistically, we can 
expect 5 of our 10 sources to be z >6.8 star-forming galaxies. 

The foregoing analysis, whilst exhaustive, is however, statistical in nature. Our ap- 
proach has been to treat all candidates as equally possible and to determine the level of 
foreground and spurious contamination as a fraction, without commenting on the nature of 
each individual source. Important information is contained in the similarity or otherwise of 
the candidate's morphology in the various NICMOS bands and the confidence with which 
we see no optical detection at the location of each candidate. For the z-band drop-outs, 
a marginal detection is permitted in the red wings of the F850LP filter, but any hint of a 
signal at the location of the candidate in shorter wavebands would give cause for concern. It 
is then a matter of judgement whether to rely primarily on the photometric redshift solution 
(§3.4.1) or to override such information and reject a candidate after visual inspection. 

Concentrating again on the z-band drop-outs, two sources, A2219-zl and A2390-zl, are 
resolved and satisfy the more rigorous color cut in Figure [3j One might therefore imagine 
these are particularly robust candidates. A2219-zl has a detection in the WFPC2 F702W 
filter very close to the location of the NICMOS image. Although reasonably significant (27.1 
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or 4.5a) it does not influence the photometric solution (Figure 6). Nonetheless, it does raise 
doubt about this candidate. 

In the case of A2390-zl, there is no detectable optical signal down to F555W= 27.6 
(2.0a) F814W= 27.0 (2.5a) and F850LP= 27.8 (2.0a) but the candidate's morphology differs 
somewhat between the NIC3 F110W and F160W filters. This may reflect the presence of 
two sources, one or both of which is at high redshift, or a genuine structural difference in 
the bands, for example result of line emission. 

The only other source worth commenting on is A2667-z2 for which there is a marginal 
WFPC2 F606W detection at the position of the NIC3 source (28.25, or 2.0a). This is 
reflected in the fact that the photometric redshift solution is fairly ambiguous for this source 
(Figure 6). 

Finally, we note as in Section 13.11 that the J — H color of most of our z drop-outs are 
significantly bluer than predicted for a normal SED at z ~7. Shifting the sources to lower 
redshift would not significantly resolve this interesting observation. 

Concerning the J-band drop-outs, neither are particularly compelling. A2667-J1 has a 
similar morphological difference between the F160W and faint F110W image and A2219-J1 
has a marginal detection (27.0, 3.5a) in the ACS F850LP filter. 



4. Spectroscopic Follow-up 

Although deep imaging with HST has delivered candidates whose photometric redshifts 
lie beyond z ~7, to date there has been not a single spectroscopic c onfirmation despite heroic 



efforts. A case in point is the z ~6.8 lensed system in Abell 2218 (IKneib et al.ll2Q04j ) which 
was the subject of 9.2ks exposure with the LRIS optical spectrograph and a marathon 33ks 
exposure with the NIRSPEC infrared spectrograph. A marginal continuum was seen but no 
e mission lines wer e detected. This contrasts with the successful detection of Lyman a in IOK- 
1 Jive et alJl20o3 ) in 31ks. The latter source has an implied star formation rate of ~ 10 M G 





yr _1 whereas the source in Abell 2218, when magnified, was expected to have an observed 
line flux equivalent to an unlensed system with a star formation rate of 2.6x25 ~ 50 M & 
yr _1 . A tantalizing explanation for the non-detection of Lyman a in the object in Abell 2218 
is preferential damping by neutral hydrogen in lower luminosity sources. 

The presence of the Lya line provides a critical feature for confirming the nature of can- 
didate high-redshift galaxies. However, the line is relatively easily attenuated and therefore 
may well be obscured in actively star forming galaxies. Therefore, the absence of the line does 
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not provide evidence that high-redshift candidates are false. Nevertheless, we conducted an 
ambitious spectroscopic campaign at Keck for some of our candidates. We naturally hoped 
that we might also secure the first spectroscopic verification of a z >7 source. 



4.1. Observations 



We used the Near InfraRed SPECtrograph fNIRSPEC. McLean et all (|l998h ) on the 10 



m Keck II Telescope to follow-up the majority of our candidates in the window 0.964-1.120 
fim, corresponding to the redshift range 6.9-8.2 for the Lyman-a line (1216 A ). Observations 
where conducted in several runs between January and September 2007 and we secured good 
data for 7 z-band drop-outs in total (see Table [5]). 

We used a 42" long and 0.76" wide slit, offering a resolving power of R ~ 1500 and used 
dithered exposures of 10 minutes each. We adjusted the dithering distance (in the range 3-8") 
in each case (third column of Table \5§ , to prevent overlap between a candidate and another 
bright source. Occasionally it was possible to observe two candidates simultaneously. For 
two candidates where we expect multiple images, the location of the expected counter-image 
(see Sect. 13.51) was used to optimize the slit positions (Figure [T2l. 



The NIRSPEC spectra were reduced following the flat-fie lding, sky-s u btracti on, dis 



tortion corrections and flux calibration procedures described by IStark et al.l (l2007bl ). These 
reduction techniques ensure an improved removal of the sky background by subtraction prior 
to resampling. We observed standard stars each night and used these to flux-calibrate the 
final spectra and determine the associated variance and hence the 5a limiting line flux. Each 
position was observed for about 3.5 to 4 hours in total, yield a limiting line flux of ~ 3 x 10~ 18 
erg cm -2 s _1 in regions of minimum OH contamination (Figure [T41) , assuming a line width 



a a ~ 300 km s 1 , as measured in the well-studied lensed system by Ellis et al.l (120011 ). For 
a different line width value, this sensitivity would vary as ^fo~^. 

We inspected each reduced 2-D spectrum for faint emission lines at the position of the 
candidate and, where relevant, that of the counter-image. No significant signal was detected 
for any of the candidates. 



4.2. Implications 

While the outcome of our spectroscopic campaign is certainly disappointing, the pres- 
ence of strong OH lines in z-band means that our limiting Ly-a flux (~3x 10~ 18 erg cm -2 
s _1 ) applies only across 50% of the observed wavelength range. Thus we would only expect 
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partial success even if all of our sources had intense emission lines. Nonetheless it is infor- 
mative to consider what the absence of any emission might mean given the star formation 
implied by our continuum detections. We can convert our flux limit into a rest-frame equiv- 
alent width using our HST photometry. For our candidates, we find a typical upper limit of 
W x ~ 5 - 20 A . 



Stanway et al.l (120071 ) have recently studied a sample of faint Lyman-a emitters at z 



6.0 selected from their photometry to be z-band drop-outs. They found a tail of high values 



for th e W\ distribution compared to similar studies undertaken at z ~ 3.0 (jShapley et al. 



20031 ). They attribute this evolution to a tendency for stronger line emission in intrinsically 
faint sources. 58% of their sample has W\ > 25 A rest-frame. Including the lower and 
upper limit measurements of W\ from this sample, which most likely contains lower redshift 
contaminants, the fraction is lowered to 34%. 

Assuming no evolution in this distribution from z ~ 6 to z ~7-8, we would expect in 
the optimistic case (58% value) ~ 4 Lyman-a emitters in our spectroscopic campaign, prior 
to considerations of the OH spectrum. The probability that all 4 objects lie in a region of 
the spectrum contaminated by OH emission is thus (0.5) 4 ~ 0.06, which is low. Even if only 
4 of the 7 candidates we examined were at high redshift (based on our statistical estimates 
given in §3), we should expect to detect the emission for ~ 3 cases. Here there would be 
only a 12 % probability of each one being occulted by OH emission. In the most pessimistic 
case from lStanway et al.l 120071 (34%), we would expect to detect only one source, with 50% 



probability of OH contamination. 

Thus, as in the case of the z ~6.8 source in Abell 2218, the absence of emission in 7 
candidates is somewhat surprising. Assuming a significant fraction are at high redshift as 
discussed in §3, this may be an important indication of the evolution in the intergalactic 
medium above z ~6. Regardless of the cause, it adds to the challenge of making progress in 
verifying high redshift candidates. 

Reconc iling the above wit h the abundance of intrinsically faint Lyman a emitters 



claimed by IStark et al.l (l2007bl ) is admittedly difficult. Should the bulk of the drop-out 
population at z >7 continue to reveal no emission, this would suggest a moderate neutral 
fraction that would challenge the transparency of the IGM at z ~8-10 implied by the pres- 
ence of feeble Lyman a emitters. The enigma simply reinforces the importance of continuing 
to attempt the detection and verification of line emission in very faint sources. 
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Discussion 



In the foregoing we have described a concerted effort to quantify the abundance of low 
luminosity star forming galaxies conducted in parallel to a simil ar spectroscopic ca mpaign 
which has examined the abundance of z ~8-10 Lyman a emitters (jStark et al.ll2007bl ). That 
study claimed that if even a small fraction of the candidates is truly at high redshift, a 
significant contribution to reionization is provided by low luminosity galaxies. In a similar 
manner, recognizing the limitations of our small samples, we now examine the luminosity 
function at z ~7.5 and the possible contribution that our lensed drop-out sources may make 
to cosmic reionization. 



5.1. Number Densities and the Source Luminosity Function 

The intrinsic area of sky (i.e. that in the source plane) covered by the NICMOS images 
is strongly dependent on the geometry of the critical lines (or caustics), which varies from 
cluster to cluster. Furthermore, multiple images occur in the central regions, duplicating the 
corresponding source plane area. The result of both effects is a smaller survey, reduced by a 
typical factor of ~10 in the source plane, with an increased depth whose value varies across 
the field of view. 

In order to derive the source density of our z-band drop-outs, and to compare our 
results with those conducted in blank fields, we used the lensing models for each cluster to 
compute the sky area effectively observed in the source plane, down to a given intrinsic AB 
magnitude. We assumed our survey covers the image plane down to the measured 5<r depth 
Jno ~ 26.2 in the central NICMOS region, and scaled according to the relative exposure time 
near the edges. We also removed ~ 10% of the NICMOS area affected by bright galaxies. 
We supposed an average redshift of z = 7.5 to compute the magnification factors. 

Errors in the magnification factors estimated from the lensing models will affect the 
source plane areas and unlensed magnitudes in an opposite way. For an individual cluster, 
the typical uncertainty is about 10% in the resulting area. This error is even smaller for the 
entire sample of six clusters. 

Our total surveyed area is a factor ~ 5 smaller than the UDF in the same magnitude 
range (AB < 27.7, Fig. [T5|) . However, the increased depth enables to reach AB~ 28 — 30 in 
this area. A very similar result is found in the case of J-band drop-outs, assuming z=9. 

We used the estimated color selection contamination factor, f cont , and the selection 
completeness factor, / comp (Sect. I3.1l and Table [3]), to correct each z drop-out individually to 
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derive intrinsic number densities. Because of the strong variations in the magnification factor 
across the NICMOS field of view, we corrected for observed completeness by computing the 
completeness factors in intrinsic (unlensed) < J + H > magnitude. This combines both the 
observed completeness factor, S comp , given in Tableland the surface reduction in the source 
plane. Error bars in the number densities were computed using Poisson noise estimates. We 
present the cumulative UV luminosity function of the z-band drop-outs in the magnitude 
range 27.0 < AB < 30.0 in Fig. ~ 



As discussed, it is likely that 5 out of our 10 sources are at high redshift. Accordingly, in 
Fig. [161 we randomly selected 100 x 5 sources from our sample to take into account object- 
to-object variations in the magnification factor and used this to estimate more realistic error 
bars. Errors in the individual magnification factors vary between 0.05 and 0.2 mag. When 
compared with the luminosity bins we used (0.5 mags), this effect has very low significance 
on the results, and the errors are mainly recovered when randomly choosing 5 candidates 
from the sample. 



F or comparison, we overplot in Fig. dn]the luminosity function found by lBouwens et al. 
( 120061 ) in the UDF, including redshift evolution betwee n z = 6.0 and z = 7. 5 assuming the 
observed size scaling as (1 + z)~ l for fixed luminosity ([Ferguson et al.ll2004f ) . Lik ewise, we 



Bouwens et al. 



(120081 ) from an 



overplot the best Schechter function fits recently claimed by 
analysis of their sample of z-band drop-outs. Not surprisingly, there is no overlap between 
these blank field measures and our, much deeper, lensed survey. All that can be said is that 
our results, which probe more than ~2.5 magnitudes fainter are marginally consistent down 
to AB~ 28.5, and higher by ~ 0.3 — 0.6 dex at fainter luminosities. 



5.2. Contribution to Cosmic Reionization 

We finally investigate whether the likely abundance of low luminosity sources found 
in our survey could make a significant contribution t o cosmic reionization . The approach 



we use is somewhat similar to the one described by IStark et al.l (l2007bl ). who estimated 
the comoving number density of sources necessary to keep the intergalactic medium (IGM) 
reionized under reasonable assumptions, and compare those to the abundances derived from 
candidate Lyman-a; emitters at high redshift. 

In our case, we can estimate the star formation rate density, measured in individual 
objects from their UV rest-frame luminosity, after applying the same completeness correc- 
tions described in Sect. 15 . 1 1 We converted the intrinsic (unlensed) J + H magnitudes of our 
candidates into a UV luminosity, L 150 o, and infer the related star formation rate (SFR) by 
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adopting the iKennicuttl (119981 ) calibration. All z-band drop-outs span the range SFR~0.1- 
1.0 M yr -1 , thus the overall star-formation rate observed yields the contribution of low 
star- formation rate sources to the entire star-formation rate density Psfr- 



We used the iMadau et al.l (119991 ) formalism to estimate the amount of star-formation 
necessary to keep the IGM reionized at a given redshift. One important factor in this 
calculation, that would modify the efficiency of star-forming sources to reionize the IGM, 
is the H n clumping factor C, defined as C =< n 2 HII > / < n HII > 2 with N HII being the 
density of ionized hydrogen. This factor measures the inhomogeneity of ionized hydrogen in 
the IGM which will likely increase between z — 10 and z = 6 due to the growth of structure. 
Assuming an IMF with a Salpeter slope with stellar masses ranging from 1 M Q to 100 M Q , 
and a solar metallicity Z=0.02 , the photon budget from star-forming sources necessary to 
reionize the IGM can be written as: 



p SFR ~ (0.031M© yr" 1 Mpc" 



esc,rel 

0.5 



-i 



C 

To 



l + z 
8.5 



where f eS c,rei is the escape fraction of ionizing photons. We assumed an escape fraction 
fesc,re/=0.5 in our calculations. Howeve r, values as low as f esCjre /=0.02 have been measured in 
z ~ 3 galaxies by IShapley et al.l (120061 ). Lower escape fractions would increase the amount 
of star formation necessary to reionize the Universe, so adopting f esc re /=0.5 gives us a lower 
limit on Psfr- On the other hand, top-heavy IMFs and differences in metallicities would 
make galaxies produce more ionizing pho tons per star-format i on ra te, but this effect is less 
significant than variations in f esCjre i and C. IBolton fc Haehnelt 



(2007 ) have critically reviewed 



possi ble values for the clumping factor C . Many authors (jBunker et al.ll2004j ; iBouwens et al. 
20061 ) have assumed C = 30 , but much lower values are predicted from radiative transfer 
simulations (llliev et al.ll2006l find C < 2 at z > 11). 



Fig. El illustrates the star formation rate densities obtain ed by integrating dow n to a 
given SFR for the two extreme luminosity functions derived by IBouwens et al.l (120081 ) when 
fitting their number densities at higher luminosities (equivalent to SFR> 1.0 M yr -1 ). We 
overplot on this figure the contributions derived from Eq. [1] with clumping factor varying 
between C = 2 and C = 30. 

As suggested before, luminous galaxies do not produce enough star formation to reionize 
the IGM at these redshifts, even when a low C = 2 clumping factor is assumed. Our 
sample enables us to test whether lower luminosity galaxies help to solve this discrepancy. 
Combining the source density from our NICMOS sur yey with those claimed at SFR> 1.0 M 
estimated by integrating the luminosity function from IBouwens et al.l (120081 ). we can compare 
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the photon budget down to SFR~ 0.1 M yr _1 with the amount of star-formation rate 
necessary to reionize the IGM. The results, shown in the case of 5 objects randomly chosen 
from our sample of 10 ^-dropouts, suggest a contribution compatible with that necessary for 
reionization for clumping factors in the range 2 < C < 10. 

6. Conclusions 

The overall goal of this project has been to constrain the abundance of low-luminosity 
star forming galaxies at z ~ 7 — 10, selected as z- and J-band drop-outs in the fields of 
6 lensing clusters observed with ACS and NICMOS onboard the Hubble Space Telescope, 
and the IRAC camera onboard the Spitzer Space Telescope. We summarize our results as 
follows: 

1. We have identified 12 high redshift candidates (10 z-band drop-outs and 2 J-band 
drop-outs) according to carefully-determined photometric selection criteria. These are 
located in 5 of the 6 clusters and span the observed magnitude range Jno —25-26. 
Each is typically magnified by 1.5 to 4 magnitudes. 

2. Based on a comprehensive set of tests, we estimate the fraction of sources that might 
represent spurious detections and the extent to which low mass stars and low redshift 
interlopers may contaminate our sample. Collectively, these tests suggest that around 
5 of of our 10 z drop-outs are possible high redshift z >7 objects. 

3. By stacking the available ACS, NICMOS, ground-based i^-band and IRAC images, we 
investigated further the averaged properties of our lensed z drop-outs. We find a UV 
spectral slope f3 ~ 2.8ig2 5 similar to that of higher luminosity candidates from the 
UDF. Such a slope suggests a very young stellar population with little reddening and 
strengthens our case that the bulk of our candidates are high redshift sources. 

4. We searched for possible counter-images for our candidates based on the most recent 
lensing models for each cluster. Unfortunately, our results are inconclusive. Many 
of the counter-images either lie outside our NICMOS field or are close to foreground 
sources. We fail to detect a counter-image in two apparently clean cases but a further 
two sources may be potentially merging on the critical line. 

5. We undertook follow-up spectroscopy with NIRSPEC for 7 of our 10 z-band drop-out 
candidates in the hope of seeing confirmatory Lyman a emission. No emission was 
found in any candidate (or its counter-image location) to a flux limit corresponding 
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to 3 x 10~ 18 ergs s _1 cm -2 in the clean part of the OH spectrum. One explanation is 
possible evolution in t he Lyman-q re s t-fram e equivalent width distribution, compared 
to previous results by lStanway et al.l (120071 ) at z ~ 6, such as might be expected if the 
neutral fraction rises with redshift. Such a deduction would be difficult to reconcile 
with t he presence of intrinsically-faint lensed Lyman a emitters at z 



10 Stark et al. 



2007b 



6. Our inferred luminosity function at z ~ 7.5, after correcting for contamination and 
incompleteness, is marginally consistent with an extrapolation of available constraints 
at brighter luminosities, with a slightly higher normalization by 0.3-0.6 dex. If even a 
modest fraction of our sources are at high redshift, our results strengthen the suggestion 
that sources with star formation rates ~ 0.1 — 1.0 M yr _1 contribute significantly to 
cosmic reionization. 



As we approach the era of JWST and the ELTs, the outcome of our project in lensed 
fields has been to provide new constraints on the faint part of the luminosity function at 
z > 7, which confirmed the trends seen at higher luminosities. Despite being restricted 
to a small field after demagnification in the source plane, we expect that these results will 
be readily confirmed by the upcoming Wide Field Camera (WFC3) on HST, combining 
extremely sensitive infrared channels with a field of view much larger than NICMOS. 
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Cluster 


RA 


Dec 


z 


Mass model 


N mutt (N z ) 


Abell 2218 


248.95625 


66.21444 


0.176 


Ehasdottir et al. (2007) 


37 (26) 


Abell 2219 


250.08541 


46.70833 


0.226 


Smith et al. (2005) 


14 (6) 


Abell 2390 


328.40086 


17.69603 


0.228 


Swinbank et al. (2006) 


11(5) 


Abell 2667 


357.91387 


-26.08541 


0.233 


Covone et al. (2006) 


10 (5) 


C10024+16 


6.65122 


17.16060 


0.390 


Kneib et al. (2003) 


9(5) 


C11358+62 


209.96069 


62.51808 


0.328 


Franx et al. (1997) 


5(3) 



Table 1: Lensing Cluster Sample: N mu i t : number of multiple images used in the lens model 
(N z number with spectroscopic redshifts) 
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Table 2. Imaging Data:. For a given cluster each entry presents the instrument, filter, 
HST/Spitzer program ID, date of observation, exposure time, and final image quality 

(depth and completeness). 



Cluster 


Filter 


Program 


Date 


Exp. time 


Depth (5a) 


Completeness 


A2218 


ACS435 


9717 


Aug04 


7048 


27. 7 '4 






ACS475 


10325 


Aug04 


5640 


27.95 






ACS555 


9717 


Aug04 


7048 


27.79 






ACS 625 


9717 


Aug04 


8386 


27.93 






ACS775 


10325 


Aug04 


9285 


27.73 






ACSsso 


9292/9452/10325 


Apr02 / Aug02 / Aug04 


19630 


27.32 






NIC3uo 


9452/10504 


Apr03/Dec05 


8446 


26.26 


25.93 




NIC3i60 


9452/10504 


Apr03/Dec05 


10559 


26.76 


26.07 




NIRC K 




Jul02 


7200 


25.5 






NIRCk 




J11IO6 


13620 


25.9 






MOIRCSk 




May07 


18000 


26.1 






IRAC3.6 


83/20439 


Dcc03/Oct05/Dec05 


37700 


25.5 






IRAC4.5 


83/20439 


Dcc03/Oct05/Dcc05 


37700 


25.5 




A2219 


WFPC2 70 2 


6488 


Aug99 


14400 


27.00 






ACSgso 


10504 


Apr06 


8374 


26.75 






NIC3no 


10504 


May06/Jun06/Jun07 


9216 


26.23 


25.97 




NIC3i6o 


10504 


May06/Jun06/Jun07 


11519 


26.73 


26.14 




NIRCk 




Jul06 


22980 


26.3 






MOIRCSk 




Aug06 


17550 


26.1 






IRAC3.6 


83 


Fcb04/Mar05 


2400 


23.9 






IRAC4.5 


83 


Fcb04/Mar05 


2400 


23.9 




A2390 


WFPC2555 


5352 


Dcc94 


8400 


26.6 






WFPC2 8 i4 


5352 


Dec94 


8400 


26.2 






ACSgso 


9292/10504 


May02/May06 


8847 


26.82 






NIC3no 


10504 


Jun06/Jul06/Jun07 


9470 


26.27 


25.89 




NIC3i60 


10504 


Jun06/Jul06/Jun07 


11839 


26.54 


26.06 




MOIRCSk 




May07 


15900 


26.0 






IRAC3.6 


83 


Jun04/Nov04 


2400 


23.9 






IRAC4.5 


83 


Jun04/Nov04 


2400 


23.9 




A2667 


WFPC2450 


8882 


OctOl 


9600 


26.26 






WFPC2 606 


8882 


OctOl 


4000 


26.94 






WFPC2 8 i4 


8882 


OctOl 


4000 


26.11 






ACSsso 


10504 


Jul06 


8765 


26.70 






NIC3no 


10504 


Aug06 


9343 


26.22 


25.93 




NIC3i6o 


10504 


Aug06 


11711 


26.51 


26.01 




IRAC3.6 


83 


Dec03 


2400 


23.9 






IRAC4.5 


83 


Dcc03 


2400 


23.9 




CL0024 


ACS435 


10325 


Nov04 


6435 


27.67 






ACS475 


10325 


Nov04 


5072 


27.81 






ACS555 


10325 


Nov04 


5072 


27.47 






ACS625 


10325 


Nov04 


8971 


27.75 






ACS775 


10325 


Nov04 


10144 


27.67 






ACSs50 


10325 


Nov04 


16328 


27.28 






NIC3no 


10504 


Aug06 


9472 


26.20 


25.9 




NIC3i6o 


10504 


Aug06 


11840 


26.60 


26.0 




MOIRCSk 




Aug06 


21600 


26.1 
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Table 2 — Continued 



Cluster Filter 


Program 


Date I 


]xp. time 


Depth (5cr) Completeness 


IRAC3.6 


64 


Dec03 


2400 


23.9 


IRAC4.5 


64 


Dec03 


3600 


24.1 


CL1358 ACS435 


9717 


Apr04/May04 


5440 


27.70 


ACS475 


9717 


Apr04/May04 


5470 


27.96 


ACS625 


9717 


Apr04/May04 


6800 


27.77 


ACS 775 


9717 


Apr04/May04 


10144 


27.49 


ACS850 


9717 


Apr04/May04 


16328 


27.13 


NIC3iio 


10504 


Dec05 


9216 


26.34 25.92 


NIC3i6o 


10504 


Dec05 


11519 


26.60 26.10 


IRAC3.6 


83 


Jan04/Jun05 


2400 


23.9 


IRAC4.5 


83 


Jan04/Jun05 


2400 


23.9 



Candidate 


R.A. 


Dec. 


-2850 (AB) 


Jiww(AB) 


H i6oiv (AB) 


K (AB) 


H(mags) 


Scomp 


fcont 


j comp 


A2218-zl 


248.9713 


+66.2071 


> 27.32 


26.1+0.13 


25.9+0.11 


> 25.7 


1.9 


0.59 


0.14 


0.66 


A2218-z2 


248.9507 


+66.2150 


> 27.32 


26.2+0.18 


26.0+0.11 


> 25.7 


2.7 


0.50 


0.19 


0.72 


A2219-zl 


250.0803 


+46.7071 


26.3 + 0.15 


24.7+0.05 


25.3+0.06 


> 25.7 


3.6 


0.84 


0.0 


0.95 


A2390-zl 


328.4130 


+17.6905 


> 26.82 


25.5+0.12 


26.1+0.12 


> 25.6 


3.5 


0.68 


0.17 


0.79 


A2390-z2 


328.4001 


+ 17.6962 


> 26.82 


25.8+0.15 


25.8+0.10 


> 25.6 


1.8 


0.65 


0.27 


0.51 


A2667-zl 


357.9119 


-26.0949 


26.7 + 0.36 


25.9+0.15 


26.1+0.18 


N/A 


1.6 


0.59 


0.34 


0.51 


A2667-z2 


357.9153 


-26.0826 


26.7 + 0.42 


25.7+0.12 


25.6+0.11 


N/A 


2.0 


0.75 


0.20 


0.68 


CL1358-zl 


209.9714 


+62.5128 


> 27.33 


26.3+0.17 


26.1+0.12 


N/A 


1.9 


0.43 


0.15 


0.63 


CL1358-z2 


209.9521 


+62.5108 


> 27.33 


26.2+0.13 


26.7+0.28 


N/A 


4.0 


0.43 


0.15 


0.72 


CL1358-z3 


209.9549 


+62.5187 


> 27.33 


26.3+0.17 


26.6+0.19 


N/A 


4.0 


0.35 


0.14 


0.66 


A2219-jl 


250.0900 


+46.7040 


> 26.7 


> 26.3 


25.0+0.05 


> 25.7 


4.0 








A2667-jl 


357.9136 


-26.0869 


> 26.7 


> 26.5 


25.1+0.08 


N/A 


3.6 









Table 3: High Redshift Candidate Photometry: total magnitudes and corresponding magni- 
fication assuming z = 7.5 for z-band drop-outs and z = 9.0 for J-band drop-outs. Each z 
drop-out entry is followed by its observed completeness, selection contamination factor and 
selection completeness. Upper limits correspond to 5 a. 



Dropouts Zg50LP JllOW Hl60VK K IRAC3.6^ m IRAC4.5 Mm 

(all 10) 28.59±0.21 25.72±0.14 25.70±0.14 

(all 4 with K-band) 28.99±0.32 25.71±0.14 25.81±0.12 > 26.2 

(all 6 with IRAC) 29.10 ±0.23 26.11±0.16 26.29±0.16 > 25.0 > 24.8 



Table 4: Stacked Photometry of the z-band drop-outs. Upper limits correspond to 5 a. 
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Run 


Candidate (s) 


Dither 


Exposure Time 


Seeing 


Notes 


Jan 2007 


CL1358zl and z2 


3" 


9 x 600s 


0.5 




May 2007 


CL1358zl and z2 


3" 


15 x 600s 


0.5-0.6" 






A2219zl 


5" 


13 x 600s 


0.5-0.8" 










12 x 600s 


0.5" 






A2218zl 


5" 


18 x 600s 


0.9" 


candidate + counter-image 








13 x 600s 


0.5" 






A2390zl 


5" 


4 x 600s 


0.8" 










12 x 600s 


0.5" 




Sep 2007 


A2390z2 


8" 


21 x 600s 


0.5" 


candidate + counter-image 




A2667zl 


6" 


17 x 600s 


0.5" 





Table 5: Log of the spectroscopic observations performed on 7 of the z band drop-outs. From 
left to right: epoch of observation, candidate name, spatial dithering between individual 
exposures, expoure time, seeing conditions. For two z drop-outs, we managed to observe the 
predicted location of the counter-image at the same time as the candidate. 



Ref. 


z 


$* (Mpc~ 3 ) 


M* (mag) 


a 


Bouwens et al. (2006) 


6 


1.4e-3 


-20.25 


-1.73 


Bouwens et al. (2008) 


7.4 


l.le-3 


-19.80 


-1.74 






1.7.8e-3 


-19.60 


-1.4 






8e-4 


-19.90 


-2.0 



Table 6: Best-fit Schechter parameters of the high redshift UV luminosity function, from 
earlier results found in the UDF, and overplotted in Fig. [T6j For each work we give the 
normali zation the absolut e magnitude M* at the exponential cutoff and the faint-end 
slope a. iBouwens et al.l (120081 ) explore three possible evolutions of the L* and $* parameters 
for different fixed slopes a. 



A AAA A AA 

0,0 " I I I I I I I I I I I I I I I I I I I I I i l l I ll 



2 3 4 5 
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Fig. 1. — Distribution of magnification factors ji (in magnitudes) for the survey, as predicted 
from the mass models assuming a point source at very high redshift [z > 7). The vertical 
axis represents the observed surface area in each 0.1 mag (0.04 dex) magnification bin. The 
peak of this distribution indicates our typical magnification factor is ~ 2.0 mags. Filled 
triangles mark the individual magnification factors for the 10 z drop-out and open triangles 
that of the J drop-out candidates. 
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Fig. 2. — Example of BCG subtraction in the clusters A2390 and A2667. Each panel is 
45"on a side. (Left to right): ACS image without subtraction, BCG-subtracted images in 
ACS, NICMOS and IRAC respectively. Subtraction residuals appear within ~ 1 arcsec of 
the center in the ACS and NICMOS case, ~ 2 arcsec for IRAC. The critical line for a high 
redshift source is show by the black curve. 
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Fig. 3. — Color-color diagram used for selecting high redshift z-band drop-out candidates. 
The final set of candidates is presented in red (se e §3.3). Color tracks represent predicte d 
colors of Hubble sequence galaxies (colored tracks) (jColeman et al.lll980l : lKinney et al.lll996l ). 
Thick tracks assume no extinction, thin lines show the effect of including a selective extinction 
of Ay = 1.0 magnitudes. T he observed location of L and T dwarfs is shown as a cross-hatched 
region (IKnapp et al.ll2004j ). Two possible z — J color selections (> 0.8 and > 1.25) are shown 
(see text for details). 




. , 1 , 1 , 1 1 1 , . 

2 4 6 8 10 

Redshift 

Fig. 4. — Optimizing the selection of high redshift z-band and J-band drop-out candidates 
using near-infrared colors (dotted lines) as compared to the expected colors of Hubble se- 
quence galaxies (as in FigJSj). 
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Fig. 5. — Testing the color-selection of z-band drop-outs. (Left) Completeness in the color- 
selection (fcomp), as a function of the J or H magnitude (assuming a flat spectrum in AB) 
and the (z — J) color-break. Values represent the fraction (%) of simulated objects whose 
photometry satisfies our color-selection criteria as indicated in the color bar beneath. (Right) 
As for the left panel but referring to the contamination fraction f CO nt- 
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Fig. 6. — NICMOS image of the cluster CL1358 (left) compared with a noise image (right) 
in order to estimate the fraction of spurious sources in our sample (see text for details). 
Except in the vicinity of the edges and central cores of bright objects (which are masked out 
by applying a simple threshold), the noise properties of the two images are very similar. 
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Fig. 7. — Snapshot images our final catalog of candidate drop-outs. (Left to right): HST 
images in the optical (ACS or WFPC2), in the ACS/F850LP (z) band, in the detection image 
(sum of F110W and F160W bands), in the NIC3/F110W and the NIC3/F160W images. To 
the right are the observed SED with the overplotted best fit template for HyperZ, and the 
redshift probability distributions (see Sect. 13.4.11) . Shown separately are the MOIRCS K 
band images.). 
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Fig. 8. — Stacked images of the z-band drop-outs. First row: ACS and NICMOS images 
of all 10 sources. Second row: stacked images for the 4 drop-outs with if-band imaging 
(rightmost image). Third row: stacked images for those 6 sources with unconfused IRAC 
data. 
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Fig. 9. — SED of the z-band drop-outs derived from the stacked photometry (see Fig. 
EJ. (Left to right): all 10 drop-outs with ACS/NICMOS data, drop-outs with useful 
ACS/NICMOS/MOIRCS data, drop-outs with useful ACS/NICMOS/IRAC data. In each 
case, the best template found with HyperZ over < z < 10 (red curve) or < z < 3 (blue 
curve) is shown. 
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Fig. 10. — Redshift probability distribution estimated using HyperZ for the stacked SEDs 
of the z-band drop-outs shown in Fig. [9j The black, red and blue curves refer respectively 
to all 10 dropouts, those four with i^-band imaging, and those 6 sources with useful IRAC 
upper limits. All three solutions are consistent with a mean population redshift of z ~ 7.35. 
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Fig. 11. — Distribution of half-light radii measured by SExtractor in the photometric 
catalogs, as a function of the total J band magnitude. The blue dashed curve corresponds 
to the measured for bright non- saturating stars. The 10 z-band drop-outs are shown as 
red ellipses. Two objects are unresolved (filled ellipses) whereas four appear resolved (open 
ellipses.). The rest cannot be reliably categorized. 
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Fig. 12. — Location of the z-band and J-band drop-outs with respect to the high redshift 
(z = 7.5) critical line for each cluster field (blue curve). Ellipses with a "CI" label mark 
the position (and estimated error) of the brightest counter-images. The adopted NIRSPEC 
follow-up slit position angles are overplotted as black rectangles. 
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Fig. 13. — A pair of magnified optical drop-outs identified in the NICMOS images. The 
theoretical location of the z = 1.8 critical line (right panel) confirms this source to be a 
low- z contaminant. 
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Fig. 14. — Determining the limiting flux from the NIRSPEC follow-up campaign. The plot 
shows the 5 a limiting emission line flux versus wavelength (and inferred Lyman a redshift) 
for a typical 3.5 hours integration. 
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Fig. 15. — Survey characteristics: the intrinsic (unlensed) surface area sampled in the source 
plane down to a given limiting magnitude for each cluster (thin red lines) and for all six 
clusters (thick black line). The upper scales give the corresponding absolute magnitudes as- 
suming z ~ 7.5. The green dashed lines illustrate the areas sampled in the abs ence of lensing. 
The blue dotted line shows the eq uivalent survey parameters for the UDF (IBouwens et al. 
2004J : bouwens fc IllingworthlhoOfil ). 
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Fig. 16. — The cumulative number density of star-forming galaxies at z ~ 7 — 8 as a function 
of the effective (unlensed) magnitude. Datapoints and error bars correspond to the range 
of densities resulting when randomly selecting 5 candidates from our sample and adopting 
Poisson errors (open diamonds, offset for clarity). In the most pessimistic case, where no 
sources are at high redshift, the implied upper li mit is shown by the th ick red curve. We 
over plot the best fit lumin osity functions found by iBouwens et al.l (120061 ) (light dotted line) 
and IBouwens et al.l (120081 ) (bold dashed lines) in the UDF (parameters in Table [6]) 
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Fig. 17. — The cumulative comoving density of star formation rate at z ^7.5 derived f or the 
two extreme (two last entries of Table [6]) luminosity functions from iBouwens et al.l (120081 ) 
with faint end slopes of a = 2.0 (purple) or a = 1.4 (green). The constraints from the 
present survey are shown as a red circle, for the average and range of densities resulting when 
randomly selecting 5 candidates from our sample (as in Fig. [T6l) . The density necessary to 
keep the IGM reionized at z — 7.5, calculated from Eq. [I] for a range of clumping factors C 
and escape fraction f escp ^eii is shown as the dashed lines. 



